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We study an underdamped oscillator with shot-noise frequency fluctuations. The oscillator spec- 
trum is determined by the interference of the susceptibilities for different eigenfrequencies. Depend- 
ing on the parameters, it has a fine structure or displays a single asymmetric peak. For nano- 
mechanical resonators with a fluctuating number of attached molecules, the spectrum is found in a 
simple analytical form. The results bear on various types of systems where the reciprocal correlation 
time of frequency fluctuations can be comparable to the typical frequency jumps. 
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Oscillators with varying frequency are studied in many 
contexts. The frequency change underlies high-resolution 
mass sensing with nano-mechanical resonators, which is 
based on the change being proportional to the mass of a 
particle attached to the resonator [J-Ql ■ It is also used in 
dynamic atomic force microscopy and in high resolution 
magnetic force microscopy 0, 0|. In quantum systems, 
the Fock states of a vibrational mode of a trapped elec- 
tron were detected from the frequency change of a non- 
linearly coupled mode 0- Recently it was proposed to 
use such change for quantum measurements of mechani- 
cal shot noise in an optomechanical system Q . 

The change of the oscillator frequency can be often 
thought of as a jump; it occurs over a much smaller time 
than the oscillator decay time or the typical inter-jump 
interval. In addition, the jumps are random. This is 
the case, e.g., where they result from attachment and de- 
tachment of molecules (or nano-particles) to a resonator 
[l|, [^, [l3l or from transitions between well-separated en- 
ergy levels in a system coupled to the oscillator 0, d, [IH . 
The oscillator dynamics is determined by the interre- 
lation between the characteristic frequency change in a 
jump A, jump rate W, and the oscillator decay rate T. Of 
utmost interest for classical and quantum measurements 
is the range where all these parameters are small com- 
pared to the oscillator eigenfrequency luq in the absence 
of jumps, W,r, A ujQ. 

In this paper we consider the susceptibility of an os- 
cillator with respect to weak resonant driving and also 
the power spectrum of the oscillator. These characteris- 
tics are advantageous as they can be directly measured 
in the experiment. We find them for an arbitrary interre- 
lation between the relevant parameters, including F, A, 
and W. The sensitivity of the spectra makes it possi- 
ble to use them for determining the parameters and the 
statistics of the underlying jump processes. 

Random frequency jumps are noise [^, [l3| , they lead 
to spectral broadening. This is an important mechanism 
of dephasing [l2|, and the analysis below immediately 
extends to systems other than the oscillator. For oscil- 



lators, the spectral broadening due to other types of fre- 
quency fluctuations has been discussed in the literature, 
see [n, 13- lH and references therein. 

The spectral broadening should be qualitatively differ- 
ent depending on the ratio A/W. For small and compar- 
atively frequent jumps, where A/W <C 1, one can think 
of the jumps as causing diffusion of the oscillator phase. 
The resulting spectral broadening should be of the order 
of the phase diffusion coefficient A'^/W [lij . 

The region A > is more complicated. One can 
imagine the oscillator as having states with different 
eigenfrequencies ujq + An separated by ^ A, where N 
enumerates the states. The interstate jump rate is ^ W. 
Such states can be thought of as different realizations 
of the oscillator and should be distinguished from the 
Fock states in each realization. Naively, one might ex- 
pect that the oscillator susceptibility xi'^) is a sum of 
independent partial susceptibilities x(-^;w) in states N. 
They would be proportional to the state populations Pn, 
and Im xi^^^) would have a form of Lorentzians cen- 
tered at + Ajv- However as shown below, this is the 
case only in the limit of large A/W. 

The inapplicability of the picture of independent par- 
tial susceptibilities can be understood by noticing that, 
in order to resolve frequencies separated by A, one should 
measure the system for time > 1/^- Therefore for 
W > A, the frequencies uJa + An cannot be resolved. As 
we show, the susceptibility can still be formally described 
as a sum of partial susceptibilities, but the latter are no 
longer independent. Rather, the complex partial sus- 
ceptibilities are coupled, their shape is strongly changed 
compared to the A/W ^ 1 limit, and the overall suscep- 
tibility can be described as a result of their interference. 
This is somewhat similar to the physics underl ying the 
paradox of the quantum harmonic oscillator [l7H20j • 

We will consider an oscillator coupled to a thermal 
bath and driven by a resonant field i^exp(— iwi)-}- c.c, 
with |w — cjol 1^0- We will assume that the coupling to 
the bath is linear in the oscillator coordinate and weak, 
and that the density of states of the bath weighted with 
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the interaction is smooth near uiq, so that the oscihator 
decay rate is essentially the same for all frequencies wq + 
Ajv,w, cf. Ref.lii 

Frequency jumps will be considered as imposed exter- 
nally, oscillator back action on the source of the jumps 
will be disregarded. This is a good approximation in 
many cases, molecule attachment-detachment and non- 
linear coupling to an off-resonance mode being examples. 
Between the jumps and in the absence of the driving and 
coupling to the bath the oscillator is described by Hamil- 
tonian Hq = h{uJo + AN)a'a which almost stepwise varies 
in time; here and a can be defined as raising and low- 
ering operators of an oscillator with frequency oj, if one 
disregards energy corrections oc (wq + Apf — uj)'^ /uj. 

The driving field does cause transitions between 
different- iV states. The oscillator response to the field 
is thus determined by the matrix elements p{N) of the 
density operator p that are diagonal with respect to N. 
Changing to the rotating frame with the canonical trans- 
formation U{t) ~ exp(— icja^ai) and using the rotating 
wave approximation, we obtain the master equation 

p{N) ^ i{Su - AN)[a^a, p{N)] + i[F'a'< + F'*a, p{N)] 
-tp{N) + Wp{N), 5u = uj-ujo, (1) 

where F' = F/{2Mhujy^^ (M is the oscillator mass). 
The operator T describes oscillator decay due to the 
coupling to the bath and has the standard form Tp = 
r(fi-|- l){a^ap — 2apa^ + pa^a) + rn{aa^ p — 2a^ pa + paa'^) , 
where fi = [exp{hu!o/kBT) — is the Planck number. 

The operator W describes transitions between states 
N with different cigcnfrequcncies, 

WpiN) ^ ^[M^(iV-r;r)p(iV-r)-M^(iV;r)p(A^)] (2) 

r 

where r enumerates the number of states over which the 
transition is made. This model describes, in particular, 
molecule attachment-detachment where molecules attach 
to a narrow region on the nano-resonator, so that the 
oscillator frequency is determined by the total number of 
attached molecules N. If molecules do not interact with 
each other, they attach/detach one by one, and 

W{N; 1) = WNo, W{N; -1) = WN, (3) 
W{N; r) = for |r| > 1; Ajv = -A^A. 

Here, A^o is the mean number of attached molecules, 
which is determined by the externally controlled molecule 
flux; for molecules of mass mmoi, A cx m„ioiWo/A/ > 0. 
The velocity jump from a mass change is small for 
nT-moi/M ^ 1, it does not cause phase accumulation in 
time in contrast to the frequency change, and can be 
disregarded; this was also checked by simulations. 

The linear response of the oscillator to the driving 
is characterized by the susceptibility X{uj) which re- 
lates the mean oscillator coordinate to the driving force, 
{q{t)) = X{uj)F exp{—iujt)+ c.c. For uj close to utQ we 



have X{ui) = (?i/2Ma;)^/^(a)/F, where the expectation 
value of operator a is given by the stationary solution 
of Eq. (P). Setting p{N) = in Eq. (UJ, multiplying 
this equation by a and taking first the trace over the 
Fock states of the oscillator for a given N (which we de- 
note by Tro below) and then the trace over A^, we obtain 
X{io) = (2Afw)~ix(w), where 

xH = ^x(A^;a.), (4) 

N 

[F - i{Slu ~ AnMN; lo) - Wx{N; to) = ^P(iV). 

Here, x(iV;w) = Troap(iV)/F', whereas P{N) = 
Tiop{N) is the stationary probability to find the oscil- 
lator in state N. From Eq. ([1]), this probability is inde- 
pendent of the driving and dissipation and is given by 
equation 

WP{N)^0, E^^(^) = l- (5) 

Equations (j4)) describes the scaled susceptibility x{^) 
as a sum of complex "partial" susceptibilities x(A^; ^) for 
each eigenfrequency state A^. These susceptibilities are 
given by a set of linear equations. They are coupled to 
each other, and since they are complex and the phase 
relations are important, one can say that they interfere, 
with xi^) determined by the result of this interference. 

Analytical expressions for xi^) can be obtained in 
the limiting cases. We start with the case where the 
frequencies in different states A^ are strongly different, 
\An - Am' \ > F, 11/ for Af 7^ A^'. Here, the partial spec- 
tra are almost independent from each other. To describe 
x{N;uj) near resonance, w « wq + Ajv, one should keep 
only diagonal in N terms in Eq. Q, which gives 

x{N;Lu)^iP{N)[^N-iiSuj-AN)]-\ (6) 
7w = F + J2,,. ^(^; ^) i\^N - An'^nI » IN, \Suj\). 

From Eq. (jB]), Im x{N;uj) has a Lorcntzian peak at fre- 
quency Wo -f A^v. The area of the peak ttP{N) is de- 
termined by the population of state A^. The halfwidth 
of the peak depends on the oscillator decay rate F 
and the total probability to switch from state A^ to other 
states, which is given by the second term in Eq. ^ for 
7Ar. The overall oscillator absorption spectrum Im xi^) 
has fine structure with peaks described by Eq. ([6]). 

In the opposite limit of small frequency change, 
|AAr|/F — !• 0, all partial spectra have the same shape, 

x{N;Lu) = tP{N){T-2Suj)-' (|Aw|«F). (7) 

The jump rate does not affect the solution in the limit 
IAatI/F — > 0. The spectrum as a whole is Lorentzian 
centered at frequency wq. This resembles the paradox of 
the quantum harmonic oscillator. The susceptibility of 
the oscillator can be presented as a superposition of cou- 
pled partial susceptibilities corresponding to transitions 
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between neighboring quantum levels. If the oscillator is 
nonlinear, the transitions occur at different frequencies 
and the corresponding spectral lines have different width, 
but for a linear oscillator, as a result of the interference, 
all partial spectra have the same shape [l^ . 

The jump rate drops out of the susceptibility also for 
arbitrary Ajv/r as long as W ^ |A7v|,r. In this case 
the term Wx{N; w) is the leading order term in Eq. 
for x{N;uj), and to first order in W'^, 

xH = i[7-«(5^-A)]-i, A = ^AArP(iV), 
7 = r + P(Af)AjvAA/x„(7V)5„(Af)(-A«)-\ (8) 

where , and Aq are the right and left eigenvectors 
and nonzero eigenvalues of matrix W (the stability of the 
oscillator stationary state implies < 0). From Eq. 
Im is again a Lorentzian peak, but now centered at 
the average frequency wq + A and with halfwidth 7 that 
exceeds T by ~ A'^/W. 

For an arbitrary relation between the oscillator param- 
eters Eqs. ([3]) can be easily solved numerically. An ex- 
plicit analytical solution can be obtained for several mod- 
els. A simple model is where frequency jumps result from 
nonlinear coupling to another oscillator (oscillator B), 
which is in thermal equilibrium; here An = NA for small 
nonlinearity (N is the energy level number of oscillator 
B). The effect of thermal transitions between neighbor- 
ing levels of oscillator B is described by Eqs. (H]), ([2]) with 
W{N; 1) = WnsiN + 1) and W^N; -1) = W{nB + l)N 
{fiB is the Planck number of oscillator B). The power 
spectrum for the corresponding frequency jumps was an- 
alyzed in Ref. [21I and the jumps as they occur in time 
were seen for an electron in a Penning trap 0- 

One of the most interesting and important for appli- 
cations is the model of molecule attachment-detachment 
Eq. ([3|). To find x{'^) it is convenient to write x(-/V, w) as 
a Fourier transform. 



x{N-uj) 



dti 



ASlu t~ 



x{N;t), xiN;0)^iPiN). 



Then Eq. (jl]) becomes a set of homogeneous differential- 
difference equations for functions x(A^;i). It can be 
solved using the discrete Laplace transform method, i.e., 
changing to x{z;t) = z^x(A^; t). This leads to a 
first-order linear partial differential equation for x{z;t). 
Its solution immediately gives x{t) — x(A"; t)^ 



m-- 



e-"+''^^««*exp 



dte'^^'xit); A/{W ~iA) 



, (9) 



Equation ([9]) gives the susceptibility of the oscillator 
with a jumping frequency as an integral of an elementary 
function. The shape of the susceptibility strongly de- 
pends on the interrelation between the frequency change 



per jump A and the jump rate W. The analysis of this 
shape can be conveniently done by rewriting Eq. ([9]) as 

x[T -i{6uj + WNoO + k{W -iA)r\ (10) 



Formally, Eq. ([TU)) for xi'^) looks like a sum of partial 
spectra. However, functions (/)fc(w) differ from the partial 
susceptibilities x(fc;a;) introduced earlier. They are close 
only in the limit of large frequency jumps, |A| ^ W, 
as seen by comparing Eqs. ([7]) and ([T0| . In this limit 
the spectrum Im x(w) is a set of well-separated equally 
spaced Lorentzian lines with halfwidths w r + {k+No)W, 
for the model Equation (|10p describes also correc- 
tions to the lineshape due to finite A/W. 

Of interest is also the limit Ao|^|^ ^ 1 where either the 
frequency shift is small, | A| ^ W, or the average number 
of attached molecules is small, Nq ^ 1. Here, the lead- 
ing order term in xi'-^) is (t'oi^)- It gives a Lorentzian 
peak of Im x(aj) centered at Slu = -(W^^/A)A^o|CPi with 
halfwidth 7 = F W^A'^ol^l^, consistent with Eq. dH) for 
A I ^ W. The frequency-jump induced broadening can 
be comparable with F even for small Ao|^|^ provided 
F < VF. For large W/\A\ we have W\^\^ oc 1/W, the 
jump-induced broadening becomes weaker with increas- 
ing W in agreement with the picture of motional nar- 
rowing [1^. The line shift Suj « —NqA is independent 
of W for PF/|A| ^ 1 and linearly increases with Nq. 
Corrections ^ -/VqICP lead to the onset of characteristic 
asymmetry of the spectrum Im x{^)- 

The evolution of the spectrum Im x{^) with varying 
parameters is illustrated in Figs. [1] and [51 Already for 
moderate frequency jumps, A/F = 3, the spectrum may 
display a well-pronounced fine structure, with the intcr- 
peak distance w A, as seen from Fig. [T] (a). With in- 
creasing jump frequency this structure is smoothed out, 
but the spectrum has a characteristic asymmetric shape. 
Fig. [1] (b). With further increase of W/ A the spectrum 
becomes narrow; the effect of frequency jumps is clearly 
seen in the shift of the peak away from ujq . 

Figure [5] demonstrates that, even for comparatively 
small A/F, not only the position of the spectral peak, 
but also its shape are sensitive to the mean number of 
attached molecules Nq. This can be helpful in the anal- 
ysis of experimental data, since A^o can be changed by 
varying the influx of molecules. 

The results on the susceptibility give also the power 
spectrum of the oscillator in the absence of resonant 
driving Q{uj) = 7r~^ Re dtexp{iujt){q{t)q{0)). Func- 
tion Q{lo) near its peak can be found from Eq. ([Ij us- 
ing the quantum regression theorem, which reduces the 
problem to equations of the form of Eq. The solu- 
tion can be obtained directly or by thinking of the fre- 
quency fluctuations as having thermal origin, and then 
Q{uj) fa (n + \){h/2TiMuj) Im x{'^) for \'-^ — ^ '^o- 
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FIG. 1. Color online. The spectrum of the oscillator with 
attaching and detaching molecules (nano-particles) for differ- 
ent scaled attachment rates W/T. The average number of 
attached molecules is A'o = 3, and A/F — 3. Data points 
show the results of numerical simulations obtained using the 
Gillespie algorithm [2^ . 
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FIG. 2. Color online. The evolution of the oscillator spec- 
trum with the varying average number of attached molecules 
A^o for a comparatively small frequency change per attached 
molecule. The parameters are A/F = 0.3, W/T = 0.1. 



The above analysis can be extended to two-level sys- 
tems with jumping inter-level spacing. If the occupation 
of the excited state can be neglected, the susceptibility is 
described by the set of linear equations (U) . The problem 
becomes more complicated if the internal dynamics of the 
system that causes frequency jumps is important. An ex- 
ample is jumps due to nonresonant coupling to a driven 
cavity mode [23[ . Here, partial susceptibilities of the sys- 
tem with jumping frequency have off-diagonal compo- 
nents with respect to the states of the system causing 
the jumps, hence instead of xi-^'i ^) oiie should consider 
x(A^, TV'; ui). If the dynamics of the jump-causing system 
is Markovian in slow time and the system of interest is in 
the ground state, one obtains from the full master equa- 
tion a set of linear equations for x(iV, N'; oj). This set is 
convenient for numerical analysis. 

The results of this paper provide an insight into the 
spectral broadening from frequency jumps and a general 
method of describing it. They show that, by studying 
coupled partial susceptibilities, one can follow the evolu- 



tion of the spectrum from well-resolved fine structure, for 
comparatively large and rare jumps, to a dephasing-type 
broadened single peak, for small and frequent jumps. The 
spectrum of a nanoresonator with attaching and detach- 
ing molecules or nano-particles is found in the explicit 
form. The susceptibility is sensitive to the parameters 
and the mechanism of the frequency jumps. 
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